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Abstract 

The pattern selection problem in binary-mixture convection in an extended chan- 
nel with a lateral through-flow is presented. The through-flow breaks left-right 
parity and changes pattern dynamics dramatically. The problem is studied based 
on computer simulation of the complete set of hydrodynamic equations (Oberbeck- 
Boussinesq approximation) in the two-dimensional rectangular channel with aspect 
ratio T = 12 and convection-suppressing lateral boundary conditions. A wide vari- 
ety of new dynamical patterns is obtained, discussed and classified. 
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Convection in a horizontal layer of binary fluid mixture heated from below is a paradigm for 
oscillatory pattern forming systems [0]]. An important generalization of the conventional type of 
the problem is the convection with a lateral through-flow. This generalization provides the unique 
opportunity to study pattern formation in a system with broken left-right parity. The symmetry 
breaking caused by the throw-flow lifts degeneracy of the problem and for this reason must affect 
the pattern dynamics dramatically. Thus, study of this type of convection sheds a light on generic 
aspects of pattern formation in systems with broken parity, which attaches especial importance to 
the problem. 

The linear stability analysis (LSA) of the spatially uniform conducting state performed in 
Ref. [Q] reveals that the through-flow splits the single Hopf bifurcation point, corresponding to 
the convection threshold without the trough-flow, into two Hopf bifurcation points for upstream 
(UTW) and downstream (DTW) traveling waves, respectively. It is important that the values of the 
thresholds as well as their relative positions with respect to each other depend on the value of the 
Reynolds number (Re) for the through-flow. Nonlinear effects are studied numerically in recent 
publication 0], where a number of interesting results is obtained. However, the small value of 
aspect ratio employed in Ref. [Q] does not allow the authors to answer the most appealing question 
"How does the trough-flow affect nonlinear stage of pattern formation and pattern selection in an 
extended channel, where the pattern structure is not imposed by the lateral boundary conditions?" 

The first attempt to answer this question is made in the present Letter. To this end we perform 
numerical integration of the standard set of two-dimensional (one horizontal coordinate x and 
vertical coordinate z, where ^-axis is antiparallel to the gravitational acceleration) hydrodynamic 
equations in the Oberbeck-Boussinesq approximation. The equations may be found, e.g., in Ref. 
[t3p and are not written here explicitly. We employ the usual non-dimensional variables, where 
length is measured in units of the layer width d and time in units of heat diffusion time d 2 / x (here 
X is the thermometric conductivity). 

The aspect ratio T is chosen equal to 12. The boundary conditions at the upper and lower 
confining surfaces are rigid, isothermal and impermeable. The lateral boundary conditions for the 
temperature and concentration fields at the inlet (x = 0) yield the basic linear profiles, while at the 
outlet (x = T) we require vanishing of the corresponding s-derivatives. As for the lateral boundary 
conditions for the velocity field, they impose the same Poiseuille profile both at the inlet and outlet. 
These boundary conditions simulate the ones, which may be imposed in real experiment, cf., e.g., 
the boundary conditions for similar problems with a through-flow employed in Refs. MM. In 
addition, suppression of the convection close to the in- and outlet allows us to minimize influence 
of the lateral boundaries on the pattern selection process. The latter has especial importance due 
to possibility of advection with the trough-flow perturbations generated by the inlet into the bulk 
of the channel. 

The Prandtl number (Pr = v/%, where v stands for the kinematic viscosity) is chosen equal 
to 10, the Lewis number (L = D/x, where D is the concentration diffusion coefficient) is 0.01. 
Such a choice of the constants corresponds to commonly used in experiment room temperature 
ethanol-water solution. For the separation ratio ip = StC (1 — C )f3/a, which is the measure 
of coupling between the temperature and the concentration gradients coursed by the Soret effect, 
the moderate value ij) = —0.1 is chosen. Here a and (3 are the thermal and solute expansion 
coefficients, St stands for the Soret coefficient and Co is the mean concentration of the solution. 
The role of control parameter plays the reduced Rayleigh number r, which at the fixed values of 
the material constants of the fluid, actually, equals AT/AT C , where AT designates the temperature 
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difference between the bottom and the top of the layer and AT C corresponds to the threshold of 
instability of the conductive state of the pure fluid without the through-flow (Re = 0, ip = 0) 
against infinitesimal perturbations. In simulations discussed in the present Letter r varies in rather 
a narrow interval close to r = 1.2. The trough-flow runs along the positive direction of x-axis 
and remains the same in all the simulations. The corresponding Reynolds number is 0.08 [|p. At 
these r and Re the characteristic values of the velocity of the shear flow are about 15-20% of the 
amplitude of the saturated convection velocity, i.e., the through-flow is always weak. 

The SIMPLE code [0] with 242x22 grid points and temporal step 5 • 1CT 4 is employed for the 
integration. The simulation gives rise to the following results. 

We begin with r = 1.20. This value of r lies close to the thresholds of instability of the 
conducting state against infinitesimal perturbations for both types of waves UTW and DTW [§] . 
The initial conditions are the conductive state with the temperature field perturbed by random 
noise with the amplitude about 10 ~ 4 . 

The simulations show that initially the conductive state becomes unstable against both the types 
of the waves. The corresponding pattern consists of superimposed, counter-propagating UTW and 
DTW. At the discussed initial stage of the instability the profiles of hydrodynamic variables for 
each type of the waves obey the relations following from the LSA (we will name these waves 
linear waves), being proportional to 

expfr 1 '' 2 '* + /jF' d x + i(k u,D x - u u ' D t)), (1) 

and similar to those in binary mixture without the trough-flow [Q]. Here the superscripts U and D 
stand for up- and downstream propagating waves and uj u < 0, [i u < 0, uj d > 0, fi D > 0. The 
wavenumbers for both types of the observed waves k u ' D w 3.2, are quite close to the values of the 
critical wavenumbers (the ones maximizing the real parts of the growth rates for UTW and DTW 
at the corresponding thresholds), given under the specified conditions by LSA [|10|]. As for u u,D , 
the expression obtained in Ref. [Q] says that in the case under consideration for both types of the 
waves 



V' D w ±6.47 + 41.9ite, (2) 



where uj^ = u u,D (k c ) at the instability threshold and signs plus and minus are related to down- 
and upstream waves respectively. The two terms in r.h.s. of Eq. @ have clear physical mean- 
ing, the first (±6.47) is associated with the Hopf bifurcation without the through-flow, while 
the second describes the drift-effect caused by the through-flow. The observed values of uj u ' d 
(uj u px —3.32, uj d ph 9.40) also are in reasonable agreement with those given by the above 
expression ph —3.11, u D ph 9.82). 

However, the real part of the growth rate for UTW is bigger than that for DTW [@]. Accord- 
ingly, soon UTW overtake DTW and begin to suppress them. During this process the amplitude 
modulations caused by superposition of the two types of waves become weaker and weaker. Fi- 
nally, DTW are suppressed completely and oscillatory growth of the amplitudes of convective 
modes is transformed into purely exponential. In other words, the system evolves to the state with 
the entire channel filled with UTW, who still obey Eq. ([I]). It happens at t « 3.5, when the maximal 
amplitude of the ^-component of the velocity field is about 10~ 2 . Thus, the first nonlinear effect 
observed in the simulation is the suppression of DTW, which agrees with the results obtained in 
Ref. [|3p for small aspect ratio. 
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Then, the growth of the amplitude of UTW brings about the next nonlinear effects, which 
change the evolution of UTW themselves. As usual [[]]], when the amplitude becomes big enough 
nonlinear stabilization occurs. The corresponding nonlinear terms affect both the real and imag- 
inary parts of the growth rate. Thus, when the amplitude reaches the vicinity of its quasi-steady 
value, determined by the balance between linear instability and nonlinear stabilization, the fre- 
quency of UTW begins to increase sharply [the frequency for the linear UTW is negative, see 
Eq. (g), so the increase of the frequency means decrease of its absolute value]. It gives rise to 
slowing down of the phase velocity of the waves and finally to their reversion, see Fig. |T]c. 

In the left half of the channel the scenario is different due to influence of the convection- 
suppressing conditions drifted from the inlet with the trough-flow. Here the exponential growth of 
the waves' amplitude suddenly, before the amplitude achieves the quasi-steady value, is changed 
with sharp decay. The convection in this half of the channel is suppressed and the fluid returns 
practically to the conductive state. The case is somewhat similar to that for convection in pure 
fluid with through-flow, where highly nonlinear states may be observed only at a certain distance 
from the inlet [0I],|120. 

Regarding the reversed waves, in the end of this stage they transform into a quasi-steady wave 
train with pronounced left boundary propagating to the outlet. However, actually these waves 
remain UTW, which are just drifted downstream due to advection with the trough-flow. The con- 
centration mode for these waves is shifted about one-quarter wavelength to the left of the w mode, 
where w stands for z-component of the velocity field in the midplane, so that counterclockwise 
circulating rolls have higher concentration level than that for the clockwise, see Fig. [j]d. Such a 
phase shift is typical to UTW, while for DTW the phase shift has the opposite sign. Another evi- 
dence of the UTW nature of these waves is their small phase velocity v p h ~ 0.43, which is smaller 
than both the maximal velocity of the trough flow = 1.20 and its average velocity Vfi = 0.8. 

At small Reynolds numbers we can write v ph (r,Re) « v ph (r, 0) + v' ph (r,0)Re, where the 
second term in the r.h.s. of this expression may be regarded as the drift velocity of the pattern 
Vdr- Generally speaking, the value of Vdr should be different from both v'J 1 ^ and Vfi due to the 
problem nonlinearity. To obtain Vd r the following approach is employed. For a given pattern we 
switch off the trough-flow (replacing the lateral boundary conditions for the velocity field with 
rigid, impermeable) and observe the dynamics of the same pattern but without the through-flow. 

As soon as the trough-flow is switched off the described quasi-steady waves change the di- 
rection of propagation and begin to travel upstream with the phase velocity v p h ~ —0.56, which 
yields Vdr the value Vd r ~ 0.99 very close to the drift velocity given by Eq. (|2|) for the linear 
waves. Bearing in mind all mentioned above, we will name these nonlinear quasi-steady waves 
drifted downstream false downstream traveling waves (FDTW). 

Since DTW are suppressed, the reversion of the left marginal wave in the UTW wave train 
creates on xt -plane a region shielded for perturbations from past to enter. The conductive state in 
this "restricted area" is destroyed with growth of small perturbations generated by the boundary 
between this region and the one filled with FDTW. It gives rise to secondary linear UTW propa- 
gating from the boundary of the FDTW wave train to the inlet, whose profiles are well described 
by Eq. ([!]), see Fig. [T]c. 

One by one the FDTW reach the outlet and vanish there. It results in contraction of the domain 
filled with FDTW and expansion of the one filled with linear UTW, respectively. It could be 
expected that finally all FDTW disappear by the outlet but, actually, it is not the case. Suddenly 
the amplitudes of FDTW begin to decay very fast. The decay yields deceleration of FDTW in 
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the laboratory coordinate frame and eventually change the direction of their propagation, i.e., the 
inverse transformation from FDTW to UTW occurs. The entire channel becomes filled with linear 
UTW and then the cycle (UTW -> UTW + FDTW -» UTW) is repeated periodically with the 
period about 24 time units, see Fig. [j]c. Possible qualitative explanation of this re-linearization 
of the convection may be as follows. The narrower the domain filled with FDTW the weaker its 
stability. When the left boundary of the wave train advances close enough to the outlet, the wave 
train cannot stand against convection-suppressing influence of the lateral boundary conditions, 
which eventually results in decrease of the waves' amplitude, i.e., to the re-linearization. 

To examine how the described nonlinear evolution depends on the Rayleigh number the simu- 
lations are repeated at different values of r. However, under the specified conditions the bifurcation 
from the conductive to convective state in binary mixture is subcritical [[[]] and the dependence of 
pattern structure on the Rayleigh number must exhibit a hysteretic loop. In this case to avoid 
jumps to another branch of the bifurcation curve the patterns obtained as the asymptotic state at 
the previous value of r are employed as the initial conditions for the new values of r. 

The simulations show that in the region from the lowest inspected value of r = 1.16 up to r = 
1.21 the pattern formation basically corresponds to the scenario discussed above. The saturated 
amplitude of the quasi-steady waves and the corresponding frequency both occur monotonically 
increasing functions of r, see Figs. [T]-|2|. The phase velocity of the quasi-steady waves vanishes at 
r « 1.17. Accordingly, FDTW do not exist at r smaller than this value, cf. Fig. [T]a and b. 

Since the effect of the sudden decay of convective modes definitely is connected with the 
convection-suppressing lateral boundary conditions, it becomes more pronounced with decrease 
of r resulting in general degradation of the convection. In accordance with that the width of the 
decay region increases, while the fraction of the period, when the quasi-steady nonlinear waves 
exist, decreases with decrease of r, cf. Fig. |T] c, b and a. Alternatively increase of r eventually 
eliminates the inverse transition FDTW — > UTW at all, see Fig. [T]e, f. 

Increase of the Rayleigh number from r = 1.21 to r = 1.213 brings about a new bifurca- 
tion. FDTW are accelerated. Both the amplitude of FDTW and their phase velocity increase and 
finally the phase velocity becomes equal to v p h ~ 1.20, which is bigger than the drift velocity, 
Fig. [l]g. Thus, FDTW transforms into true DTW. The transformation is accompanied with the 
corresponding change of the phase shift between the concentration and w modes, so that now 
counterclockwise circulating rolls have smaller concentration level than that for the clockwise, see 
Fig. [j]h. In the same time linear UTW in the left half of the channel do not undergo the sudden 
decay any more. Eventually the system evolves to a steady state with left part of the channel 
filled with the linear UTW, which obey Eq. ([I]), and right part with the nonlinear steady DTW. The 
boundary between these two domains remains practically immobile. 

Another important feature of the problem revealed in simulations is strict wavenumber and 
frequency selection. Every value of the Rayleigh number yields a certain unique pair of k and u for 
each of the discussed types of the waves. Any change of r triggers transient process, which results 
in the corresponding changes of k and to. The dependence of k and w on r for nonlinear quasi- 
steady UTW, FDTW and DTW is presented in Fig. ^. It is important to stress that the imposed 
boundary conditions prevent any phase pinning or wavenumber discretization, so the observed 
wavenumber selection is the generic, intrinsic feature of the problem, which is not related to any 
external factor. 

Note that transition UTW — > FDTW does not change the smoothness of curves k(r) and u>(r), 
while transition FDTW — > DTW does, see Fig. ^. An interesting feature of the phenomenon is 
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that the functions k(r) and u>(r) for UTW and FDTW are much sharper than those for DTW, cf. 
the behavior of these function blow and beyond the transition point. At decrease of r a narrow 
hysteretic loop for inverse transition DTW — > FDTW is observed. DTW remain stable until r = 
2.09. Further decrease of the Rayleigh number to r = 2.08 triggers transformation DTW — > 
FDTW with formation of a pattern similar to that shown in Fig. [j]c. 

The last issue to be discussed is the global stability of the observed waves. According to the 
results of Ref. [^p for an unbounded in the lateral dimension channel steady nonlinear UTW should 
be unstable at employed in our simulations values of ip = — 0.1, Re = 0.08 and any value of r. The 
stable patterns observed in our simulations do not contradict to this result because their stability 
may be associated with the size-effect at the moderate value of T = 12 used in the simulations. 
To elucidate this question simulations with bigger value of T are required. These simulations are 
under way and will be reported elsewhere. 
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FIG. 1. Spatiotemporal plots for z-component of the velocity field in the midplane w at different values 
of the Rayleigh number and profiles of w (black), concentration c = (C — Co)/?/ (ctAT) (red) and temper- 
ature 8 = (T — To)/ AT (blue) fields in the midplane at some characteristic moments of time, a, b, c, e — 
asymptotic, periodically repeated in time patterns at r = 1.16 (a), r = 1.18 (b), r = 1.20 (c) and r = 1.21 
(e), respectively, d — profiles of w, c and 9 for spatiotemporal plot c at two moments of time indicated in c 
as white lines, the upper and lower panels on d correspond to the upper and lower lines on c. f — analogous 
profiles for e. g — transition FDTW — > DTW initiated by abrupt raise of r from 1.21 to 1.213 at t = 0. 
h — profiles of w, c and 6 for plot g. The profiles in the lower panel (the lower white line on g) correspond 
to FDTW, while those in the upper panel (the upper white line on g) already to DTW. For more details see 
the text. 
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FIG. 2. The dependence of the wavenumber k (o) and oscillation frequency uj (o) on the reduced 
Rayleigh number r for nonlinear quasi-steady UTW, FDTW and DTW. Note the discontinuity at r = 1.213 
caused by transition FDTW — > DTW. Closed symbols designate the boundary of a hysteretic loop observed 
during decrease of r. Full straight line in the plot corresponds to to = 0, dashed to the drift contribution of 
the through-flow for the linear waves [second term in the r.h.s. of Eq. (Q)]. 
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